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Abstract 

This study explores how the metal materials of the applicator influence the dose distribution when performing 
brachytherapy for cervical cancer. A pinpoint ionization chamber, IVlonte Carlo code MCNPX, and treatment planning 
system are used to evaluate the dose distribution for a single lr-192 source positioned in the tandem and ovoid. For dose 
distribution in water with the presence of the tandem, differences among measurement, IVICNPX calculation and treatment 
planning system results are <5%. For dose distribution in water with the presence of the ovoid, the IVICNPX result agrees 
with the measurement. But the doses calculated from treatment planning system are overestimated by up to a factor of 4. 
This is due to the shielding effect of the metal materials in the applicator not being considered in the treatment planning 
system. This result suggests that the treatment planning system should take into account corrections for the metal materials 
of the applicator in order to improve the accuracy of the radiation dose delivered. 
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Introduction 

Statistics from the World Health Organization indicates that 
cervical cancer is the second most common cancer among women 
after breast cancer [1]. In addition to invasive surgical procedures, 
the main treatments for cervical cancer now include chemother- 
apy, external-beam radiotherapy, and brachytherapy [2,3]. 

The 5-year survival rates of patients with Stage IB or IIA 
cervical cancer who receive external-beam radiotherapy together 
with brachytherapy are about 80~90%, which is as high as for 
patients treated with hysterectomy together with pelvic lymph- 
node dissection [4]. Moreover, radiation therapy is better than 
surgical treatment for controlling the tumor in patients with later 
stages of cervical cancer [5] . Brachytherapy delivers a radiation 
dose to the tumor over a short distance. The steep dose gradient 
result in the tumor receiving a high dose within a short time and 
reduces the doses received by adjacent healthy tissues. Brachy- 
therapy doses are currently calculated based on the TG-43 report 
of The American Association of Physicists in Medicine [6] . The 
parameters and dose calculation formulas provided in that report 
are based on the assumption of a pure water environment. 
However, the change in the brachytherapy dose distribution in 
actual treatments due to the inhomogeneous media of the human 
body is an important issue clinically [7] . 

An applicator is used to position the radioactive source within 
the patient's body when performing brachytherapy treatment for 
cervical cancer. Usually the applicator has a metal material 



covering its outside and might even also contain a metal structure 
within the ovoid. Henschke applicators have been widely used 
inpatients with cervical cancer in recent years [8-1 1], but the dose 
perturbation caused by such apphcators is not currently evaluated 
in treatment planning system. The European Society for 
Therapeutic Radiology and Oncology(ESTRO) suggested consid- 
ering the effects of tissue inhomogeneity and source shielding 
during treatment planning system [12]. Watanabe used thermo- 
luminescent dosimeter and MCNP4Acode to perform dose 
measurements and dose calculations for a Henschke applicator, 
with the results revealing uncertainty between the MCNP4A- 
calculated value and the actual dose measured [13]. 

The present study aims to determine how the dose distribution 
is influenced by shielding materials contained in the ovoid foranlr- 
192 source inside the Henschke applicator. Monte-Carlo N- 
Particle eXtended (MCNPX) code used in this study to simulate 
the dose distribution of the Henschke applicator, and a small 
ionization chamber is used for measurements. The measured 
values are compared with the results of MCNPX and the results of 
the treatment planning system. The deviation of dose delivery 
caused by dosimeter displacement is also explored. 

Materials and Methods 

lr-192 source and Henschke applicator 

The lr-192 source used in this study is 3.5 mm long with a 
diameter of 0.6 mm. It is covered by stainless steel with an outer 
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Figure 1. Structural diagram of tlie lr-192 source and its stainless steel outer cover. 

doi:1 0.1 371 /journal.pone.01 04831 .gOOl 



diameter of 1 . 1 mm. The source is attached to a remote after- 
loading machine (MicroSelectron HDR, Nucletron, Netherlands) 
by a stainless steel wire. Fig. 1 shows the geometric diagram of the 
source. Fig. 2 shows the Henschke applicator, which comprises of 
a tandem and ovoid son its two sides. The applicator functions as 
the translation pathway for the source during cervical cancer 
treatment. The ovoid contains tungsten alloy as a shielding 
material with a density of 17.0 g/cvc?, the composition of which is 
91% tungsten, 4.5% nickel, and 4.5% iron. Kodak X-Omat V 
films are used to capture the shielding structure and layout inside 
the ovoid of the Henschke applicator. MCNPX code is used to 
simulate the Henschke applicator and Ir-192 source. With the aid 
of MCNPX Visual Editor, the applicator model set up is shown in 
Fig. 2. 

Experimental setup of dose distribution measurement 
using Henschke applicator 

In this study a water phantom (MT-IOO-T, MED-TEC, USA) 
with dimensions of 30.5 cmx38 cmx38 cm is used to perform 
measurements with Ir-192 located at a single position to determine 
how the metal materials of the Henschke applicator influences the 
dose distribution. This experiment used a 0.01-ccpinpoint 
ionization chamber (CCOl, IBA Dosimetry, Germany) to measure 
the dose distribution for the ovoid and tandem. The X, Y, and Z 
axes correspond to left-right, front-back, and up-down directions, 
respectively. The measurement area is 0 mm in the X direction, 
from 15 mm to 45 mm in the Y direction, and from —40 mm to 
40 mm in the Z direction. 



Dose distribution calculation using MCNPX 

In this study the differences among the ionization-chamber- 
measured, treatment-planning, and MCNPX-calculated values are 
compared. The reliability of the measured values can be 
confirmed by MCNPX results. 

Comparison of the MCNPX F6 and *F8 tallies 

MCNPX was developed by the Los Alamos National Labora- 
tory (LANL) [14]. Two tally methods F6 & *F8 are used to 
calculate the MCNPX photon dose for comparison. F6 tally is 
based on the assumption that electronic balances exists in the tally 
region, and that an electron would lose its energy instead of 
vmdergoing electron transport tat the position where a photon and 
electron collide. F8 tally considers the energy deposition of an 
electron on its migration path when a photon and electron collide, 
rather than the instantaneous energy deposition at the position of 
the collision. The F8 tally can reflect more truly the actual 
migration of radiation particles and provides more accurate 
calculation results, but the computing time is longer. 

Both the F6 and F8 tallies are used to calculate the Ir- 1 92 dose 
distribution in the water phantom to see the dose differences 
between these two tally methods. The simulation includes only the 
Ir-192 source, and not the Henschke applicator. The Ir-192 
photon spectrum was obtained from Brookhaven National 
Laboratory [15]. Tally spheres with diameter of 1 mm are placed 
at distances of 15 mm and 25 mm on the Y axis. Tally spheres 
with diameter of 3 mm are placed at distances of 35 mm and 
45 mm on the Y axis. For number of particles equals 5x10, the 
statistical errors of less than 3%. In addition, the F6 tally is used to 
calculate anisotropy functions and radial dose functions for 
comparison with the results of Williamson and Li [16]. 



3,60 




Figure 2. Structural model of Henschke applicator for Ir-192 source for MCNP calculation, visualized by MCNPX Visual Editor. 

doi:1 0.1 371/journal.pone.01 04831. g002 
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Figure 3. Diagram of the setup for simulating dose variation caused by imprecision of measurement-position MCNPX tally spheres 
are placed at different distances from the source. Positive/negative signs indicate positions away from or close to the source, respectively. 
doi:1 0.1 371/journal.pone.01 04831. g003 



(a)i.ioo 



(C) 



1.000 



sr 0.800 



;|::::t::i:::p:iii::l 



♦ Y-axis = 15 mm | 
■ Y-axis = 25 mm i 
A Y-axis = 35 mm 

• Y-axis = 45 mm 



-10 10 
Z-axis (mm) 




(d) 




— MCNPX, r = 1.0 cm (Water) 
- MCNPX, r = 2.0 cm (Water) 
■ - MCNPX, r = 3.0 cm (Water) 
■--MCNPX, r = 5.0cm (Water) 
A Williamson and Li, r = 1.0 cm 

• Williamson and U, r = 2.0 cm 
■ Williamson and U, r = 3.0 cm 

• Williamson and U, r = 5.0 cm \ 



SO 100 
6 (degree) 



120 



140 



160 



4.00 6.00 8.00 10.00 
Radial distance r (cm) 




MCNPX Y-axis = 15 mm 
«-MCNPXV.a]iis = 25 mm 

MCNPX Y-axis = 35 mm 
-•-MCNPX Y-axis = 45 

TPS Y-axis - 15 mm 

• TPS Y-axis = 25 mm 

• TPS Y-axis = 35 mm 

• TPS Y-axis = 45 mm 
I.e. Y-axis = 15 mm 

□ I.e. Y-axts = 25 mm 
A I.e. Y-axis = 35 mm 
c I.e. Y-axis = 45 mm 



Figure 4. Monte Carlos simulation, (a) Ratios of F6 and F8 tally dose of MCNPX calculation for the presence of lr-192 in the water phantom, (b) 
Comparison of MCNPX-calculated and previously reported results for the lr-192 radial dose function, (c) Comparison of MCNPX-calculated and 
previously reported results for the lr-192anisotropy function, (d) Comparison of dose distribution in water with the presence of the tandem for the 
ionization-chamber-measured value (I.C.), treatment-planning system value (TPS), and MCNPX-calculated value. 
doi:1 0.1 371/journal.pone.01 04831. g004 
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Figure 5. Relative dose differences in water witli the presence of the tandem, (a) IVICNPX-calculatecJ value versus treatment-planning value, 
(b) MCNPX-calculated value versus ionization-chamber-measured value, and (c) treatment-planning value versus ionization-chamber-measured value. 
doi:1 0.1 371 /journal.pone.01 04831 .g005 
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Figure 6. The ratio of dose witKi the presence of the ovoid and with the presence of the tandem. 

doi:1 0.1 371 /journal.pone.01 04831 .g006 



Monte Carlo calculation for the dose distribution in water 
with presence of the Henschke applicator 

For the dose distribution in water witli presence of tlie 
Henschke applicator, MCNPX code is used to calculate the dose 
to clarify the reliability of the ionization-chamber-measured 
values. The conditions of the simulation are made as close as 
possible to those in actual measurements. The F6 tally is used for 
dose calculations in the simulation in order to reduce the 
calculation time and statistical errors. The MCNPX-calculated, 
ionization-chamber-measured, and treatment-planning calculated 
values are compared to assess how the metal materials influence 
the dose distribution. 

Dose variation caused by imprecision of measurement 
position 

Brachytherapy is associated with steep dose gradients. To 
evaluate dose variation caused by imprecision of positioning 
during the experiment, MCNPX is applied to simulate the 
situation of an Ir-192 source at a frxed position. The tally spheres 
are at 15 mm, 25 mm, 35 mm, and 45 mm away from the source 
as shown in Fig. 3. The imprecision of measurement position are 




X-axis 



Figure 7. Photograph showing the internal structure of the 
ovoid of the Hensclii<e applicator. 

doi:1 0.1 371/journal.pone.01 04831. g007 



±0.25 mm, ±0.50 mm, ±1.0 mm, and ±2.0 mm along the Y 
axis. 

Results and Discussion 

Differences between F6 and *F8 tallies in MCNP 
calculation 

The ratio of doses calculated by using F6 and F8 tallies on the 
in MCNPX calculation is shown in Fig. 4(a). The value Hes 
between 0.97 to 1.03. That is the differences between these two 
types of tally are small, <3%. However, the computing time of *F8 
tally is ~30 times longer than that of the F6 tally. Therefore, 
performing an Ir-192 dose calculation using MCNPX, F6 tally can 
be used when the environment media is homogeneous to save 
computing time. 

Radial dose function and anisotropy function 

Figure 4(b) shows the radial dose function, i.e. the dose 
reduction profile caused by photon absorption and scattering 
while the photons travel along the transverse direction of the 
source calculated by using MCNPX. The results agree well with 
those of Williamson and Li [16]. The differences are <1.8%, less 
than the 2a statistical error (=3%) of the MCNPX calculation. 
The number of particle used in MCNPX calculation is 5x10^'. 

The anisotropy function accounts for the anisotropy of dose 
distribution around the source, including the effects of absorption 
and scatter in the medium. The calculation results are compared 
with those of Williamson and Li [16] as shown in Fig. 4(c). For 
region of 1 cm radius the differences are within 4.6% for 9<5°, 
within 2.5% for 5°<9<180°. For region outside 1 cm radius, the 
differences are within 2%. All are within 2cr of MCNPX 
calculations. 

Dose distribution in water with the presence of the 
tandem 

To evaluate the accuracy of the measured values under the 
situation with presence of tandem, this study compared the 
measurement results with the calculated results of MCNPX and 
the treatment-planning system. The measured and the calculated 
results of dose distribution in water with the presence of the 
tandem are shown in Fig. 4(d). The measurement and the 
calculated results of treatment-planning and MCNPX are all 
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normalized to dose at Y = 1 5 mm and Z — 0 mm. The MCNPX 
results shows that in Y-direction, the dose at 15 mm, 25 mm, 
35 mm and 45 mm drop from 1 to 0.36, 0.18 and 0.11, almost 
follows the inverse square law. The average dose reduction is 6.4% 
per mm in Y direction. For Y = 15 mm, the dose along z-direction 
changes from 1 at z = 0 to 0.68 at z = 10 mm. The average dose 
reduction in Z-direction is 3.2% per mm. The dose gradient is 
steeper along the Y direction than along Z direction. Hence the 
dose errors due to imprecision location would be greater along the 
Y direction. 

The comparison of this three relative dose distribution are 
shown in Fig. 5. After normalized, the dose differences between 
the MCNPX-calculated and treatment-planning results are < 
1.5"'o, within 2a of the Monte Carlo calculation. The difference 
decreases as distance from source increases. The dose diflFerences 
between the ionization-chamber-measured and MCNPX-calcu- 
lated values are <5%; the dose differences between the measured 
and the treatment-planning system values are also <5%. 

Dose distribution in water witli tine presence of the ovoid 

Current treatment planning system for brachytherapy does not 
consider the dose attenuation due to the metal inside the ovoid. 
The shielding effect of metal inside the ovoid is analyzed for both 
measured and MCNPX calculation results. The ratios of dose 
calculated by MCNPX with the presence of ovoid to dose with 
tandem present only are shown in Fig. 6. For Y= 15 mm, the 
ratio changes from 0.94 at outside region to 0.22 at inside region 
(<-l — 10 mm, as shown in Fig. 7). The effect is less severe as the 
distance from the source axis increases. The difiTerence of dose 
ratio between the measured and the MCNPX calculated values 
are ~5-8%. Therefore, if the shielding effect of the ovoid is 
neglected, the dose predicted by treatment planning system can be 
overestimated by up to a factor of ~4. 

The bladder and rectum are located within the adjacent areas 
above and beneath the Y axis. It is verified by both clinical 
measurements and MCNPX calculations that tlu- Ilcnschke 
applicator does provide shielding effect and can reduce the doses 
to the bladder and rectum. The current treatment-planning system 
which neglects this metal-shielded area on the other hand, will 
overestimate the dose of the shielded area by a factor of as large as 
4. The above result also implies that if the tumor is located behind 
the ovoid, the tumor dose will be overestimated if the shielding 
effect of the metal is not taken into consideration. The magnitude 
of the overestimation of dose depends on the dwelling time of the 

References 

1. Allrn D, Narayan K (200.^) Managing advanced-stage ccndcal cancer. Best Pract 
Res Clin Olistet Gynaecol 19: 591-609. 

2. Macdonald DM, Lin LL, Bichl K, Mutic S, Nantz R, et al. (2008) Combined 
intensity-modulated radiation therapy and brachytherapy in the treatment of 
cervical cancer. IntJ Radiat Oncol Biol Phys 71: 618—624. 

3. Tanderup K, Georg D, Pdtter R, Kirisits C, Grau C, et al. Adaptive 
management of cervical cancer radiotherapy; 2010. Semin Radiat Oncol, pp. 
121-129. 

4. Kim HJ, Kim S, Ha SW, Wu H-G (2008) Are doses to ICRU reference points 
valuable for predicting late rectal and bladder morbidity after definitive 
radiotherapy in uterine cervbc caneeri' Tumori 94: 327-332. 

5. Undurraga M, Loubeyre P. Dubuisson J-B, Schneider D, Petignat P (2010) 
Early-stage cervical cancer: is surgery better than radiotherapy? Expert Rev 
Anticancer Ther 10: 451-460. 

6. Nath R, Anderson LL, Luxton G, Weaver KA, Williamson JF, et al. (1995) 
Dosimetry of interstitial brachytherapy sources: recommendations of the AAPM 
Radiation Therapy Committee Task Group No. 43. Med Phys 22: 209-234. 

7. Hsu S-M, Wu C-H, Lee J-H, Hsieh Y-J, Yu C-Y, et al. (2012) A study on the 
dose distributions in various materials from an Ir-192 HDR brachytherapy 
source. PLoS One 7: e44528. 

8. Gifford KA, Hortonjr JL, Pelloski CE, Jhingran A, Mourtada F, et al. (2005) A 
three-dimensional computed tomography— assisted Monte Carlo evaluation of 



source in the region shielded by the metal of ovoid from the source 
location. 

Dose deviation caused by the imprecise positioning 

The dose deviation caused by imprecise mcasun^ment-position 
calculated by MCNPX is shown in Table 1. The dose deviation 
for imprecise positioning of ±0.25 mm is small, ~3%, comparable 
to the 2cTstatistical error. For imprecise positioning of ±0.5 mm, 
±1 mm, and ±2 mm, the dose deviation are 6~7%, 12~15% 
and 22~33% respectively. 

The dose deviation caused by imprecise positioning tended to 
decrease with distance from the source. Special attention should 
therefore be paid to precise positioning when measuring doses in 
areas near the source location. This observation also implies that if 
the source location relative to tumor is imprecise, dose error up to 
17% per mm would occur. Therefore, the applicator setup errors 
should be reduced to within 0.25 mm to improve the accuracy of 
the dose delivery during the brachytherapy. 

Conclusions 

This study used a pinpoint ionization chamber to perform dose 
measurements for the Henschke applicator containing a source at 
a single point. In addiation to the treatment planning system, 
MCNPX is used to simulate the experimental condition to confirm 
the measurements. For dose calculation with presence of tandem, 
the three results, after normalized, agree well with one another. 

For dose calculation with presence of ovoid, the MCNPX- 
Cedculated and ionization-chamber-measured values agree well < 
5%, the difference is within 2a of MCNPX calculation. 

The treatment planning system, which currently does not take 
into account the shielding material in ovoid, overestimation of the 
dose by up to a factor of 4 is possible. This results show that the 
shielding matericd inside the ovoid of the Henschke applicator can 
reduce the dose delivered to adjacent healthy organs. However, 
the tumor located behind the ovoid, its dose will also be 
overestimated by the treatment planning calculation. Corrections 
should be made for the presence of the Henschke applicator to 
ensure the accuracy of the doses delivered. 
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